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The paper is concerned with the geometrically non-linear theory of 6-parame- 
tric elastic shells with drilling degrees of freedom. This theory establishes a 
general model for shells, which is characterized by two independent kinematic 
fields: the translation vector and the rotation tensor. Thus, the kinematical 
structure of 6-parameter shells is identical to that of Cosserat shells. We show 
^ . the existence of global minimizers for the geometrically non-linear 2D equations 

of elastic shells. The proof of the existence theorem is based on the direct methods 
of the calculus of variations using essentially the convexity of the energy in the 
strain and curvature measures. Since our result is valid for general anisotropic 
shells, we analyze separately the particular cases of isotropic shells, orthotropic 
shells, and composite shells. 
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1 Introduction 



In recent years there has been a revived interest in 2D shell models because of uncon- 
ventional materials and extremely small aspect to thickness ratio, such as for instance 
thin polymeric films or biological membranes. For classical engineering materials and 
for non-extreme aspect to thickness ratios, available 3D FEM-Codes may readily be 
used such that the need for a truly 2D shell model does not arise anymore. However, for 
ultrathin specimens the application of a 3D constitutive law is not clear at all. In these 
extreme cases one is led to employ a 2D shell model. This paper is concerned with one 
such model, the geometrically non-linear resultant theory of shells. We consider the 
6-parameter model of shells which involves two independent kinematic fields: the trans- 
lation vector field and the rotation tensor field (in total 6 independent scalar kinematic 
variables). This theory of shells is one of the most general, and it is also very effective 
in the treatment of complex shell problems, as can be seen from the works [T5ll25l[5T| . 
among others. The resultant 6-parameter theory of shells was originally proposed by 
Reissner [55] and it has been considerably developed subsequently. An account of these 
developments and main achievements have been presented in the books of Libai and 
Simmonds [32] and Chroscielewski, Makowski and Pietraszkiewicz [14J. In this ap- 
proach, the 2D equilibrium equations and static boundary conditions of the shell are 
derived exactly by direct through-the-thickness integration of the stresses in the 3D bal- 
ance laws of linear and angular momentum. The kinematic fields are then constructed 
on the 2D level using the integral identity of the virtual work principle. Following this 
procedure, the 2D model is expressed in terms of stress resultants and work-averaged 
deformation fields defined on the shell base surface. It is interesting that the kine- 
matical structure of 6-parameter shells (involving the translation vector and rotation 
tensor) is identical to the kinematical structure of Cosserat shells (defined as material 
surfaces endowed with a triad of rigid directors describing the orientation of points). 
From this point of view, the 6-parameter theory of shells is related to the shell model 
proposed initially by the Cosserat brothers [20] and developed by many authors, such 
as Zhilin [65], Zubov [66], Altenbach and Zhilin [5], Eremeyev and Zubov [26], Birsan 
and Altenbach [TD]. Using the so-called derivation approach, Neff [3HII12] has estab- 
lished independently a Cosserat-type model for initially planar shells (plates) which 
is very similar to the 6-parameter resultant shell model. A comparison between these 
two models has been presented in the paper [11], in the case of plates. 

On the other hand, we should mention that the kinematic structure of the 6- 
parameter shell model is different from the kinematic structure of the so-called Cosserat 
surfaces, which are defined as material surfaces with one or more deformable directors 
attached to every point, see [ll[6l[71[35l[56l[53 . For instance, the kinematics of Cosserat 



surfaces with one deformable director is characterized also by 6 degrees of freedom (3 
for the position of material points and 3 for the orientation and stretch of the mate- 
rial line element through the thickness), which differ essentially from the 6 degrees of 
freedom in the 6-parameter resultant shell model. 

The topic of existence of solutions for the 2D equations of linear and non-linear 
elastic shells has been treated in many works. The results that can be found in the 
literature refer to various types of shell models and they employ different techniques, 
see e.g. [Il[2l|8l|9l[211[28l[30l[3ll[58H6l]. The method of formal asymptotic expansions is 
one method of investigation which allows for the derivation and justification of plate 
and shell models. The existence theory for linear or nonlinear shells is presented in 
details in the books of Ciarlet [T614T8] . together with many historical remarks and 
bibliographic references. Another fruitful approach to the existence theory of 2D plate 
and shell models (obtained as limit cases of 3D models) is the F-convergence analysis 
of thin structures, see e.g. [1111121 SSI 113 SB]- Concerning the geometrically non-hnear 
6-parameter theory of elastic shells, there is no existence theorem published in the 
literature yet, as far as we are aware of. In the case of linear micropolar shells, the 
existence of weak solutions has been recently proved in [22]. Existence results for 
the related (very similar) Cosserat-type model of initially planar shells have been 
established by Neff [38lll2]. In [HlllHlllS] the linearized version of this model has been 
analyzed and compared with the classical membrane and bending plate models given 
by the Reissner-Mindlin or Kirchhoff-Love theories. 

In the present work, we prove the existence of minimizers for the minimization 
problem of the total potential energy associated to the deformation of geometrically 
non-linear 6-parameter elastic shells. We emphasize that our work is not concerned 
with the derivation of the 2D shell model, but it presents existence results for the well- 
established 2D theory of 6-parameter elastic shells. It should be mentioned from the 
beginning that this model refers to shells made of a simple (classical) elastic material, 
not a generalized (Cosserat or micropolar) continuum. However, the rotation tensor 
field appears naturally in this theory, in the course of the exact through-the-thickness 
reduction of the 3D formulation of the problem to the 2D one [Mll2^ l32] . Thus, in spite 
of the above mentioned similarity with the kinematics of Cosserat shells, the material 
of the shell in the resultant 6-parameter model is described as a simple continuum 
(without any specific microstructure or material length scale). On the other hand, in 
the case of dimensional reduction of the 3D equations of micropolar shell-like bodies 
one can obtain the same 6-parameter theory with modified 2D constitutive equations, 
see e.g. [^ for the linear case and [38lll2l|67] for the nonlinear case, or one can obtain 
more complex theories as in [27]. 

For the proof of existence, we employ the direct methods of the calculus of varia- 



tions and extend the techniques presented in [3S1I12] to the case of general shells (with 
non- vanishing curvature in the reference configuration). In Section 2 we present briefly 
the kinematics of general 6-parameter shells and the equations of equilibrium. In Sec- 
tion 3 we give some alternative formulas for the strain tensor and curvature tensor, 
which are written in the direct tensorial notation as well as in the component (ma- 
trix) notation. These expressions are needed subsequently in the proof of our main 
result. In Section 4 we formulate the two-field minimization problem for general elas- 
tic shells, corresponding to mixed-type boundary conditions. Under the assumptions 
of convexity and coercivity of the quadratic strain energy function (physically linear 
material response), we prove the existence of minimizers over a large set of admissible 
pairs. Thus, the minimizing solution pair is of class i?^(a;,R'^) for the translation vec- 
tor and H^{uj, 5*0(3)) for the rotation tensor. The existence result is valid for general 
anisotropic elastic shells having arbitrary geometry of the reference configuration. Sec- 
tion 5 includes some applications of the existence theorem and discussions of special 
cases. We present a convenient way to choose the initial directors and the parametriza- 
tion of the reference surface. Then, we consider separately the cases of isotropic shells, 
orthotropic shells, and composite layered shells and we present the respective forms 
of the strain energy densities. Applying the theorem stated previously, we establish 
the conditions on the constitutive coefficients that ensure the existence of minimizers 
in each situation. This analysis shows the usefulness of our theoretical result in the 
treatment of practical problems for elastic shells. 

2 General 6-parameter resultant shells 

Consider a general 6-parameter shell and denote with the base surface of the shell 
in the reference (initial) configuration and with S the base surface in the deformed 
configuration. Let O be a fixed point in the Euclidean space and {ei, 62, 63} the fixed 
orthonormal vector basis. The reference configuration is represented by the position 
vector (relative to the point O) of the base surface plus the structure tensor Q^. 
The structure tensor is a second order proper orthogonal tensor which can be described 
by an orthonormal triad of directors {d^, cZ^, d^} attached to every point [THElj. Thus 
the reference (initial) configuration is characterized by the functions 



where thus (xi, X2) are material curvilinear coordinates on the surface . Throughout 
the paper Latin indexes take the values {1,2,3}, while Greek indexes ... 

the values {1,2}. The usual Einstein summation convention over repeated indexes 



y° : w C R2 ^ M^ 
Q° : a; C ^ S0{3) 



yO = yO(xi,X2), 



(1) 



Q\y') 




Figure 1: The base surface of the shell in the initial configuration, the base surface 5* 
in the deformed configuration, and the fictitious planar reference configuration uj. The 
orthonormal triads of vectors {cj} , {<i°} and {dj} are related through the relations 
di = Q^dPi = Rsi and = Q^e, , where is the elastic rotation field, is the 
initial rotation, and R is the total rotation field. 



is employed. We assume that the curvilinear coordinates (xi,X2) G u range over a 
bounded open domain u (with Lipschitz boundary du) of the OxiX2 plane, see Figure 
[U Let us denote the partial derivative with respect to by daf = , for any 
function /. We designate by {ai, the (covariant) base vectors in the tangent plane 
of and by nP the unit normal to given by 

aa = OaV = t: — , = . (2) 

OXa \\CLl X a2\\ 

The reciprocal (contravariant) basis {a^, a^} of the tangent plane is defined by a'^-ap = 
5^ (the Kronecker symbol). We also use the notations 

ag = = n°, a^is = ■ , a°'^ = a" ■ a'^, a = ^Jdet{a^^^^)^ > 0. 

For the deformed configuration of the shell, let y{xi,X2) denote the position vector 
(relative to O) and {di{xi, X2)} the orthonormal triad of directors attached to the point 
with initial curvilinear coordinates (xi,X2). The deformed configuration is completely 
characterized by the functions 

y = x{y'), Q'^ = di®d°G50(3), (3) 



/ 
/ 

Figure 2: Multiplicative decomposition of the total deformation gradient F = F^F^ into 
the elastic shell deformation gradient F^ and the initial deformation gradient F^ = P. In- 
terpretation in terms of reconstructed three-dimensional quantities. The elastic response is 
governed by F'^. The curved initial configuration corresponds to the intermediate stress free 
configuration in multiplicative plasticity. 

where x : S"*^ — )■ represents the deformation of the base surface and the proper 
orthogonal tensor field is the (effective) elastic rotation. The displacement vector 
is defined as usual hy u = y — . 

We mention that the role of the triads of directors {d^} and {di] is to determine 
the structure tensor and the rotation tensor of the shell, respectively. Thus, the 
directors do not describe here any microstructure of the material. According to the 
derivation procedure of the 6-parameter shell model, the kinematical fields y and 
are uniquely defined as the work-conjugate averages of 3D deformation distribution 
over the shell thickness, whose virtual values enter the virtual work principle of the 
shell (see [I1E3]). 

In view of ([1]) and (|3]), the deformed configuration can alternatively be characterized 
by the functions 

y = y{xi,X2) = x(y°(xi,a;2)), R{xi,X2) = Q'^Q^ = di{xi,X2) ® Si e 5*0(3), 

where the vector y and the orthogonal tensor R are fields defined over u. The or- 
thogonal tensor field Q'^ represents the elastic rotation tensor between the reference 
and deformed configurations [S21IS3]- The tensor is the initial rotation field, while 
R = Q^Q^ describes the total rotation from the fictitious planar reference configura- 
tion uj (endowed with the triad {ej}) to the deformed configuration S. We mention 
that the tensors Q° and R are also called the structure tensors of the reference and 
deformed configurations, respectively [111121]. The following relations hold 

Q' = RQ''', <i° = Q°ei, d, = Q'd^ = Re,. (4) 




Usually, the initial directors d- are chosen such that dg = nP and belong to the 
tangent plane of (see Remark [TU]) . This assumption is not necessary in general and 
we do not use it in the proof of our existence result. 

Let F = GradgT/ = daV ® a,"' denote the (total) shell deformation gradient tensor. 
The strong form of the equations of equilibrium for 6-parameter shells can be written 
in the form [21] 

Div, N + f = 0, Div, M + ax\{NF'^ - FN'^) + c = 0, (5) 

where N and M are the internal surface stress resultant and stress couple tensors of 
the 1** Piola-Kirchhoff type, while / and c are the external surface resultant force 
and couple vectors applied to points of S, but measured per unit area of . The 
operators Grad^ and Div^ are the surface gradient and surface divergence, respectively, 
defined intrinsically in [221131] • The superscript {■)'^ denotes the transpose and axl( ■) 
represents the axial vector of a skew-symmetric tensor. 

Let u be the external unit normal vector to the boundary curve dS^ lying in the 
tangent plane. We consider boundary conditions of the type [23|l5T] 

Nv = n*, Mh> = m* along dS^ , 

y = y*, R = R* along dS^, ^ ' 

where dS^ = a5jU(95° is a disjoint partition of 5° {dS^^ndS^ = 0) with \ength{dS^) > 
0. Here, n* and m* are the external boundary resultant force and couple vectors 
respectively, applied along the deformed boundary dS, but measured per unit length 
of dS'j C dS^ . We denote by duf and dud the subsets of the boundary curve du which 
correspond to dSj and dS^ respectively, through the mapping . 

The weak form associated to these local balance equations for shells has been pre- 
sented in [Tl|23l[32]. 

3 Elastic shell strain and curvature measures 

According to [lllEl], the elastic shell strain tensor E'^ in the material representation 
is given by 

= Q^'^Grad, y - Grad, y' = {Q^'^d^y - d^y') ® a" , (7) 

since Grad^y = d^y ® a," . It is useful to write the strain tensor E*^ in the alternative 
form 

E^ = {Q^''^day - ««) ® a" = (Q^'^S^y ® a" + n° ® a^) - (a^ ® a') 
= {Q^'^daV ® e« + nO ® eg) (e^ ® a') - 



or equivalently, since (cj (g) a*) = (a^ (S) ej) ^, 

= {Q^'^do^y ® + n° ® eg) (a^ ® e^) - II3 = C7' - Us , 
with = (Q'^'^Say (g) Cq, + n° (g) 63) (ttj (g) Bj) \ 



(8) 



where II3 = Cj is the identity tensor and U represents the non-symmetric elastic 
shell stretch tensor, which can be seen as the 2D analog of the 3D non-symmetric 
Biot-type stretch tensor |13] or the first Cosserat deformation tensor [20j page 123, eq. 
(43)] for the shelL Let us denote by P the tensor defined by 



P = ei = daV^ (g) + n° (g) 63 . 

Then, from (|8]) and (Q we get 

E'=Tj"-lis = Q"'^{day ® + Q^nO ® 63) - II3 , 



(9) 



U = Q'^^{d^y ® e„ + Q'n'' ® 63) P 



-1 



(10) 



In the sequel, it is useful to write the elastic shell strain tensors in component form, 
relative to the fixed tensor basis {ej ® ej}. Let = {Eij)^^r^ be the matrix of 
components for the tensor E'^ = EfjCi (8 ej . In general, we decompose any second 
order tensor T in the form T = TijBi (8 ej and denote by T = (^2^)3^3 the matrix of 
components. Also, for any vector v = ViCi we designate by f = ("^^1)3^^ the column 
matrix of components. 



Remark 1 The matrix of components P = {Pij)^y^r^ for the tensor defined in (|9 
be specified in terms of its 3 column vectors as follows 



P 



d2y' 



n 



3x3 



n 



3x3 



ai 



02 



n 



3x3 



can 



(11) 



We mention that the tensor P introduced in 
(deformation) gradient 

P = V0(xi,X2,X3) 



can be seen as a three-dimensional 



\X3=0 



with 



0(xi, X2, X3) = y^{xi, X2) + X3 n°(xi, Xa) 



(12) 



and it satisfies det P = a/ det Oq,^ = a > 0, so that the inverse P ^ exists. The mapping 



e-.ujxi-ll) 



has been introduced previously in [T^lT^ lS ^ and employed for 



the geometrical description of 3D shells (h denotes the thickness of the shell). 

By virtue of ([HD and ([9]), we obtain the following matrix form for the strain tensor E^ 



Q'^^d,y 



(13) 



where II3 = {Sij)ny^r, is the unit matrix. The matrix E"^ can be written equivalently as 



d2y 



E^ 



ge^O j p-1 

Q^'^F^- II3 



II3, or 



(14) 



with 



:= (^diy 
F := (^diy d2y 
F^ ■=P=[ diy^ 



(^diy d2y Q^n°^ P-^ = Q^'^ (^Vy Q^n^^ P~^ 



Q^n^^ P~^ = (^Vy Q"n°)(Ve)-^ 



d2y' 
F 



Vy 



(15) 



n 



Here, F*^ represents the (reconstructed) elastic shell deformation gradient, F is the (re- 
constructed) total deformation gradient and F" = P is the initial deformation gradient 
and we recognize an analog of the classical multiplicative decomposition F = F^F^ of 
the deformation gradient into elastic and initial ( "plastic" ) parts [371110] , see Figure [2l 
The elastic material response is defined in terms of the elastic part of the deformation, 
e.g. E^ = Q^'^F^ - II3 , cf. ([HD. 

Remark 2 Although the resultant shell model is truly a 2D theory, we may always 
consider reconstructed three-dimensional quantities. In this sense, similar to the context 
of RemarkUl the tensor F = daV^ea + Q^n^^e^, which has the components matrix 
F and which represents the reconstructed total deformation gradient, can be interpreted 
as a three-dimensional deformation gradient 



Vcp(xi,X2,X3) 



U3=0 



with 



(p{xi,X2, X3) = y{xi,X2) + xs Q''{xi,X2)'nP{xi, X2) 

= y{xi, X2) + X3 Q''{xi,X2) V @{xi,X2, 0)63. 



(16) 



Here, the mapping ip : uj x (— f,!) — i-K^ is a reconstructed 3D deformation of the 
shell, in terms of the given 2D quantities y{xi,X2) and Q^{xi,X2)- 

In terms of the total rotation R and the initial rotation Q*^, the elastic shell strain 
tensor is expressed by 

E'' = Q° {R^d^y - ® a" . (17) 

Then, we have 

= Q°[(i?^5„y-Q°''^9,yO)®eJ (e, ® a') = Q'[{R^d^y-Q'^^d^y^)®e^]p-' 



which can be written in matrix form as follows 



= Q^H P-^ with ^ := ( R^diy - Q^'^diy^ 



) . (18) 

3x3 



On the other hand, the elastic shell curvature tensor in the material description 
is defined by [HISl] 

= [Q^'^ax\{d^RR^) - axl(9„Q°Q°'^)] ® a" . (19) 



In order to write in a form more convenient to us, we use relations of the type 

Q^ax\{d^QQ^) = axl(Q^9„Q), axl(QAQ^) = Qaxl(A), (20) 

which hold true for any orthogonal tensor Q G 50(3) and any skew-symmetric tensor 
A G SO (3) (see e.g., p^). Using f pUj) in f lT^ we can write the elastic curvature tensor 
in the equivalent forms 



(21) 



or. 



= [axliQ^'^d^Q") ® e„] (a, ® e^) ^ = [axl(Q"'^9„Q^) ® e^]P~\ 
Then, the matrix of components = (^K^j^^^^ is given by 



K''= SixliQ^'^diQ^) axl(Q"'^92Q 



p-\ 



(22) 



If we express K"^ in terms of the total rotation R and the initial rotation Q^, we 



get 

= Q° [ax\{R^dc,R) - axl(Q°'^9„g°)] ® a" . 
This relation can be written as 

= K- with K := Q^qx\{R^ d^^R) ® a" 



(23) 



(24) 



a" 



Q°axl(Q°'^9„Q°) ® a" = axl(9„Q° Q°'^) 

where the tensor li' is the total curvature tensor, while is the initial curvature (or 
structure curvature tensor of 5*°). In view of (|23|1 and (El]), the matrix = {K^j) is 
given by 

= Q^LP-^ = K-K° with 



L :-- 



axl{R'^diR) - axl(Q0'^9ig°) ax\{R^d2R) - ax^Q^.r^^Q 







3x3 



= Q0( d.x\{R^diR) ax\{R'^d2R) 



p-\ 



(25) 



)p-^ 



i^o = Qo( axKQO.T^^QO) I ax^Q^.T^^gc 

In what follows, we shall use the expressions f llSp and f l25p of the elastic shell strain 
measures E'^ and written with tensor components in the basis {cj Cj}. 



Remark 3 As expected, the case of zero strain and bending measures corresponds to 
a rigid body mode of the shell. Indeed, if E'^ = and = 0, then from ([7]) and 
we obtain 

daV = Q"^ daV^ and daQ"" = 0. 
Hence, it follows that is constant and 

y = Q^lP + <^ = constant), 

which means that the shell undergoes a rigid body motion with constant translation c 
and constant rotation Q^. 

Remark 4 In the case when the base surface 5*° of the initial configuration of the shell 
is planar we may assume that coincides with u. In this situation we have ai = Bi , 
P = II3 , and the above strain and curvature measures coincide with those defined for 
the Cosserat model of planar-shells introduced in JS3i\T^- 

Remark 5 In view of (12T]) or ( 122|) . the elastic shell curvature tensor is an analog 
of the second Cosserat deformation tensor in the 3D theory, see the original Cosserats 
book [20, page 123, eq. (44)]. 



4 Variational formulation for elastic shells 

Let us denote the strain energy density of the elastic shell hj W = W{E^, K^). Ac- 
cording to the hyperelasticity assumption, the internal surface stress resultant N and 
stress couple tensor M are expressed by the constitutive equations in the form 

dW dW 

In this paper we assume that the strain energy density is a quadratic function of 
its arguments E'^ and K"^. Thus, the considered model is physically linear and geo- 
metrically non-linear. The explicit form of the strain energy function W is presented 
in plll32] for isotropic, hemitropic or orthotropic elastic shells. In general, the coeffi- 
cients of the strain energy function W depend on the structure curvature tensor 
see [21] . In [13] , the case of composite (layered) shells is investigated and the expression 
of the energy density is established. These special cases will be discussed in Section [51 
Consider the usual Lebesgue spaces [Lp{u), \\ ■ \\lp{ui)), P > 1, and Sobolev space 
[H^{uj), II ■ We denote by L^{uj,R^) (respectively i/^(a;,]R^)) the space of all 

vector fields v = ViCi such that v-i G L^ioo) (respectively Vi G H^{u)). Similarly, we 
denote the sets H\uj,R^''^) = {T = T.jSi | T^^- G H\co)} , H\co,SO{3)) = 



{T G H\u,R^''^)\T e S0{3)} , XP(a;,M3x3) = {T = T.^e, ® ej\Tij e , 
^0(3)) = {T G I ^ 50(3)}. The norm of a tensor T is defined by 

\\Tf = ti{TT'') = TijTij . 

Concerning the boundary- value problem ([5]), ([6]), we assume the existence of a 
function A{y, R) representing the potential of external surface loads /, c, and boundary 
loads n*, m* (cf. [23]). 

We consider the following two-field minimization problem associated to the defor- 
mation of elastic shells: find the pair {y, R) in the admissible set A which realizes the 
minimum of the functional 

Iiy,R)= [ W{E',K')dS-Aiy,R) for {y,R)eA, (27) 

where dS is the area element of the surface . The admissible set A is defined by 

A= {iy,R) e H\uj,W') X H\uj,S0i3)) \ y\gso = y*, R\^^, = R*}, (28) 

where the boundary conditions are to be understood in the sense of traces. The tensors 
E'^ and K'^ are expressed in terms of {y, R) through the relations ([T7|) and f l23|) . If we 
write W = W{E'', K") = W(yy, R, VR), then referring the integral to the (fictitious 
reference) domain u, the change of variable formula clearly gives 

W{E^, K^) dS = [ W{E', K') a(xi, xa) da^idxa 
so ^ (29) 

= W[Vy{xi,X2), R{xi,X2),VR{xi,X2)) det (V0(xi, X2, 0)) dxidx2, 

where a = ^/det{a^ is the notation introduced previously. The variational princi- 
ple for the total energy of elastic shells with respect to the functional fl27p has been 
presented in [23], Sect. 2. We decompose the loading potential A{y,R) additively as 
follows 

A{y, R) = Aso{y,R) + AQso{y,R), 



Asoiy,R)= [ f-udS + UsoiR), AesoXy,R)=[ n*-udl + U 

Jso ' JdS" 



(30) 



ISO ' JdS"j ' 

where u = y — y^ is the displacement vector and dl is the element of length along 
the boundary curve dS^^ . In fl30|) . Aso{y,R) is the potential of the external surface 
loads f,c, while Agso{y,Q^) is the potential of the external boundary loads n*,m*. 
The expression of the load potential functions II^o , Ilg^o : L'^{co, 5*0(3)) — )► M are not 
given exphcitly, but they are assumed to be continuous and bounded operators. Of 
course, the integrals over 5*° and dS'j appearing in (!30|) can be transformed like in (!29|) 
into integrals over u and duj , respectively. 

We mention that one can consider more general cases of external loads in the 
definition of the loading potential ( 15U]) . such as for example tracking loads. 



4.1 Main result: Existence of minimizers 

This theorem states the existence of minimizers to the minimization problem f l27l) - fl30|) . 



Theorem 6 Assume that the external loads satisfy the conditions 

f G L\uj, M^), n* G L^duf, M^), (31) 

and the boundary data satisfy the conditions 

y* e H\uj,M.^), R* e H\u,S0{3)). (32) 

Assume that the following conditions concerning the initial configuration are fulfilled: 
7/° : a; C — )■ M'^ a continuous injective mapping and 

y"^ e H\cu,m^), Q"^ e H\uj,S0{3)), (33) 

aa = d^y^ G X°°(a;, M^) ( i.e. G L°°(a;, ^ 

(34) 

det {aap{xi,X2)) > cto > ' 

where oq is a constant. The strain energy density W{E^, K^) is assumed to he a 
quadratic convex function of {E'^^K^) and W is coercive, in the sense that there exists 
a constant Co > with 

W{E^,K^) > Co{\\E^\\^ + \\K^\\^). (35) 
Then, the minimization problem f l27|) -f l30|) admits at least one minimizing solution pair 

Remark 7 The hypotheses can be written equivalently in terms of the tensor 
P = V0(xi,X2,O) as 

P G X°°(w,M3x3)^ detP>ao>0, (36) 

in view of the relations (|9]) and ( [TT]) . Since y^ represents the position vector of the 
reference base surface (which is bounded), the conditions ( l33l) i and ( l34ll i can be 
written together in the form y° G VF"'^'°°(a;, M^). 

Proof. We employ the direct methods of the calculus of variations. We show first that 
there exists a constant C > such that 

|A(j/,P)| < + V(z/,P)gA (37) 



Indeed, since G L°°{u,M.^) it follows that a = ^ydet{a^ G L°°{u). We also have 
ll-Rp = tT{RR^) = 3, V-R G 5*0(3). Taking into account the hypotheses (^T]) and the 
boundedness of 1150 and Ilg^o , we deduce from fl30|) that 

I A(t/, R) \ < \ Aso{y, R)\ + \ Aos^.{y, R)\ < Ci \\y - y°|U2(^) + C2 \\y - y'^Wmdu^f) 

+ 1 Uso{R) I + 1 nacjo(i?) I < C3||?/|U2(,) + C4||2/||hi(^) + C5 , 

for some positive constants Ck > 0. Then, the inequality (^71) holds. 

In what follows, we employ the component form of the elastic strain tensors and 
K^, written as matrices E'^ and K'^ in f|T8l) and fl25l) . respectively. Let us show next 
that there exists a positive constant Aq > such that 

ll^^'ll = ll^'ll > Ao||i/|| , \\K'\\ = \\K^\\ > Xo\\L\\, (38) 

where the matrices H = (ifjj)^^^ and L = (Ljj)^^^ are introduced in f|T8l) and fl25l) . 
Indeed, since = Q^HP-^ and Q° G 50(3) we have 

II f = II Q'HP-' f = II i/p-i f = tr[/7P-i(/7P-i)^] = tr [//(P'^P)"!//^] . (39) 
From flTT]) we deduce that 
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Inserting (HOl) into (!39l) we obtain 

II E' ||2 = a'^^H.^H,^ = a'^^H^^H^p + a'^^H^Ap + a^^^H^^H.^ , with = (f^^.Og.g , 

(41) 

since Hi^ = according to (ITS!) . In virtue of ([31]), it follows that the matrix ((3.0^)2x2 
and its inverse matrix (fl"^)2x2 ~ ('^a/s) ^ satisfy 

(a,^) G L°°(w,M'^') and (a"^) G X°°(a;, M^^'). 

Then, the smallest eigenvalue of the positive definite symmetric matrix (a°'^(xi, 3:2)) 2x2 
is greater than a positive constant Aq > and consequently 

a°'%xi,X2) VaVi3 > Agf^f^, V (xi, 3:2) G w, Vfi,f2G]R. (42) 

Using inequality (l42l) for each individual sum in the right-hand side of (l^Ti) we deduce 
that II-E'^IP > XlHiaHia = Aq ||i?|P, i.e. the inequality is proved. The proof of 



the inequality ( 135]) 9 is identical. In view of and ( ITS]) we have 

0=1 

= A^^dl R'dM'- 2(i?^9„y,Q°'W) + II Q°'Wir) 
«=1 

= a£(|| a„2/f - 2(i?^a„y ,Q°'W) + II 5«/ir), 



0=1 



where {S,T) = tr[ST'^] is the scalar product of two matrices S,T. Integrating over u 
and using the Cauchy-Schwarz inequality we obtain 

2 

II Whi.) > E (II ^"^ll'^H - 2|l do.y\\LHu.)\\ d.y' IIl^h + || dMU..)), 

a=l 

or 

II Whiu.) > (II dMlhi.) + II 92y\\U.)) -C4y \\h^^.) + , (43) 

for some positive constants Ci > 0, C2 > 0. Let us show that the functional I{y, R) is 
bounded from below over the admissible set A. By virtue of ( l29l) . ( I34l) 9 and ( 1371) we 
can write 

liy, R)>cJ II f a dxidx2 - A(2/, i?) > ^ ao|| - C ( \\y\\m(^) + l) 

and using ( H3|) we deduce that there exist the constants C^> and such that 
I{y,R) > CoaoXl (|| aiy||i.(,) + || d2y\\l.^^))-Cs\\y\\HH.)-C,, V(y,K) G ^ (44) 



with oo specified by (IMj) . We observe that the vector field y — y* ^ 
satisfies y — y* = on dua ■ Applying the Poincare-inequality we infer the existence 



of a constant Cp > such that 



my - 2/*)|li.H + my - I/*)||i.H >cp\\y-y* ||^.(,) . (45) 

Using inequalities of the type \\day\\l2^^) > {Wd^iy - ?/*)||l2h - ||5a?/*||L2H)^ and 
[5l) we find that 



ll^ii/lli^H + W^^yWhn ^ 11^ - y*\\miu) 

-2||t/ - rllz/^Mdl^irlU^H + l|52rl|L^H) + (II 5irlli2(,) + II d^yXn.))- 



From the last inequality and f l44j) follows that there exist some constants C5 > and 
Cq with 

/(y, i^) > Co aoXl Cp \\y - y*\\l.^^^ - C^y - y*\\mH + 0^, V {y, R) e A. (46) 



Since the constant Cq ao-^o Cp > 0, the function I{y, R) is bounded from below over A. 
Hence, there exists an infimizing sequence \^{y^, Rn)}^_-^ C A such that 

hm I{y^, Rn) = inf {l{y, R) \ {y, R)eA}. (47) 

In view of the conditions fl32|) we have I{y*,R*) < oo. The infimizing sequence 
\^{y^, Rn)}^_-^ can be chosen such that 

Iiy^,R^) < Iiy*,R*) <oo, Vn > 1. (48) 

Taking into account fH^ and fHS]) we see that the sequence {yn}^_i is bounded in 
H\u,R^). Then, we can extract a subsequence of {yn}^-i (not relabeled) which 
converges weakly in H'^{uj, M^) and moreover, according to Rellich's selection principle, 
it converges strongly in L\u,R^), i.e. there exists an element y G H^{u, M^) such that 

y^^y in H\uj,R^), and Vn ^ V in L^iuj,R^). (49) 

For any n G N, let us denote by and li'^ the strain measures corresponding to the 
fields [y^, Rn), defined by the relations (HZD and (ES}. We have E"^, G L2(w,M3x3) 
and let i?^, i^^ be the matrices of components in the basis {ej (g) Cj}, given by f lTSj) and 
in terms From ([27]), (ElDa , dSS]), (EZ]) and (|1H]) we get 

^^0^0 II Kl ||i,(^) < I W{El Kl) a(xi, x^) dx^dx, < I{y*, R*) + C { \\yJmH + l)- 

Since {yn}^_^ is bounded in H^{u,R^), it follows from the last inequalities that 
{K^}^^^ is bounded in L^{uj, M^xS). in view of , we deduce that {axl{RldaRn)}^ 
is bounded in L'^{oo,R^), or equivalently {daRn}^^^ is bounded in L^(a;,M^^^), for 
a = 1,2. Since i?„ G 5*0(3) we have H-RnP = 3 and thus we can infer that the 
sequence {-R„}^^ is bounded in H^{u,R^^'^^). Hence, there exists a subsequence of 
[Rn}^_^ (not relabeled) and an element R G H^{uj,R^^^) with 

Rn^R in Jf^(w,M^^^), and R^ ^ R in ^^(a;, M^xs^,^ ^gg^ 
We can show for the limit that R G 5*0(3). Indeed, since Rn G 5*0(3) we have 

II RnR^ — I13||l2(w) = II Rn{R^ — R^)\\l'^{ui) = || -R — -Rn||L2(w) > 0, 

i.e. RnR^ — > Is in L'^{u,R^^^). On the other hand, we can write 

II RnR^ — RR^\\l^{ui) = II {Rn — R)R^\\l^{u) < 3|| i?„ — R\\l2i^^^ \\ i?||/,2(^) )■ 0, 

which means that RnR^ RR^ in i"'^(a;, M^^^). Consequently, we find RR^ = ^3 
so that i?G fi"^(w, ^0(3)). 



By virtue of the relations -R„) G A and (H^ . ( ISUj) . we derive that y = y* on 
dS'^ and i? = R* on in the sense of traces. Hence, we obtain that the hmit pair 
satisfies {y, R) G A. 

Let us construct the elements E'^,K'^ G 2y^(u;,M^^^) defined in terms of the fields 
{y,R) by the relations f[T7|) and f l2^ . Then, the matrices of components E^^K^ are 
expressed in terms of the components (y, -R) by f lTS]) and fl2^ . i.e. 

= Q0( - QO'^^iyO I RTQ^y - QO.T^^y" I ) P-\ 

(51) 

ir^ = QO(axl(i?^9i/?)-axl(gO'^9igO) | axl(i?^92i?) -axl(g0'^92g°) | 0)p-i. 

Next, we want to show that there exist some subsequences (not relabeled) of {-E^} 
and {K^n} such that 

El^E^ in 2v2(a;,M3x3^)^ K^^ in ^^(u;,^^^^). (52) 

As shown above, the sequence {y^j^j^ is bounded in i/^(a;,]R^). It follows that 
{c^a?/n}^i bounded in L'^{cu, M^) and then the sequence {R^dayn}^_i bounded in 
2v^(a;,M^), since J2„ G 5*0(3). Consequently, there exists a subsequence (not relabeled) 
and an element ^„ G L'^{uj, M'^) such that 

i?l5a?/n ^ la ^ ^'(a;,^^). (53) 

On the other hand, let (f) G C[5°(a;,M'^) be an arbitrary test function. Then, using the 
properties of the scalar product we deduce 

/ [R^dayn - R^d^y) ■ ^dxidx2 

= / lf{do,yn - 9o,y) ■ (I>dxidx2 + / [R^ - R^)dayn ■ (I>dxidx2 

< \\Rn-R\\L^uj)\\dayn®(f>\\L^u)+ [day d^y) • Rcf) dxidx2 , 

since the relations (H9ll . ( I50ll and i?<^ G -L^(a;,R^) hold, and ||9q,2/,„ (g) 4>\\ is bounded. 
Thus, we get 

! [Rld^yn] ■ ct>dxrdx2 I [R^day] ■ c^dx.dx^, V0 g C^{u,R'). (54) 



By comparison of (1531) and flMl) we find €q, = R dalJ , which means that R^day^ 
Wday in X^(a;,M^) , or equivalently 

[Rldayn-Q'^'^daV') ^ {ieday-Rlday'') in 2.2(0;, R^). (55) 



Taking into account (fTSi) . (15Ti) i and the hypotheses (155]) . we obtain from (1531) that 
E" in 1v2(w,M3x3), i.e. the relation holds. 
To prove fl52l) 9 we start from the fact that the sequence {R^daRn}'^_i is bounded in 
L'^{co,M!^^^), as we proved previously. Then, there exists a subsequence (not relabeled) 
and an element G L^(a;,M'^^^) such that 

RldaRn^Ca "1 (w, M=^^=^) . (56) 

On the other hand, for any test function # G C^(c<;, R'^^^) we can write 

/ {R^daRn - R^daR, ^) dxidx2 = / {R^[daRn - d^R) , * ) dxidx2 

J U} J U] 

+ / {[Rl - R^)daRn , * > dxidxa < / (d^Rn - d^R, R$ ) dxidx2 

+ ||-Rn — -R||l2(<^) \\daRn^'^ Wl'^iuj) ^ 0; 

since R$ G X^(w,]R^^^), \\daRn^^\\ is bounded, and relations (150|) hold. Conse- 
quently, we have 



/ {R^daRn , dxidx2 — > / {R^daR, dxidx2, V * G C, 



and by comparison with ( l56l) we deduce that = R daR, i.e. the convergence 
Rld^Rn R^daR holds in L2(w,M3x3)^ follows that 

[ax\{R^daRn) - axl(i2;^9aQ°)] ^ [axl(^^9a^) - axl(i2^9«Q°)] in X2(w,M^^^), 

and from (!25|) . (!33l) . and (!5T|) 9 we derive that the convergence (152!) 9 holds true. 

In the last step of the proof we use the convexity of the strain energy density W. 
In view of fl52l) . we have 



W{E^, K^)a{xi, X2)dxidx2 < liminf / W{E^^, K^J a{xi,X2) dxidx2. (57) 

since is convex in [E'^, K^). Taking into account the hypotheses flHT]) . the continuity 
of the load potential functions Ilgo , Ilg^o , and the convergence relations (H9!l 9 and 
(!50|) 9 , we deduce 

K{y,R)= \im K{y^,Rr,). (58) 
From ([27D, ([29D, ([57D and ([58]) we get 

I{y,R) < liminf /(?;„, i2„). (59) 

n— ^oo 

Finally, the relations fHTI) and (159!) show that 

/(y,^) = inf {/(2/,il)|(t/,il)G^}. 



Since {ij, R) G A, we conclude that {ij, R) is a minimizing solution pair of our mini- 
mization problem. The proof is complete. ■ 



Remark 8 The solution fields satisfy the following regularity conditions 

if G H\uj,M.^), Re L'^{io,SO{3))r\H\io,SO{3)). 

Thus, the position vector y and the total rotation field R may fail to he continuous, 
according to the limit case of Sobolev embedding. 

Remark 9 We observe that the boundary conditions imposed on the orthogonal tensor 
R can be relaxed in the definition of the admissible set A. Thus, one can prove the 
existence of minimizers for the minimization problem (1271) over the following larger 
admissible set 

A= {{y,R) e H\u,R') X H\u,S0i3)) \ y^^^^ = y*}. 

This assertion can be proved in the same way as the Theorem 0. For a discussion of 
possible alternative boundary conditions for the field R on dud we refer to the works 

5 Applications of the theorem and discussions 

In this section we present some important special cases for the choice of the energy 
density W where the Theorem [6] can be successfully applied to show the existence of 
minimizers. 

Let us discuss first on the choice of the 3 initial directors {d^} in the reference 
configuration, i.e. the specification of the proper orthogonal tensor Q° = d° (g) . One 
judicious choice for the tensor is the following 

g° = polar(P) = polar(Ve(a;i,X2,0)), (60) 

where P = Cj = (9a!/° ® Cq, + n° (g) 63 has been introduced previously in ([9]) and 
polar(T) denotes the orthogonal tensor given by the polar decomposition of any tensor 
T. 

Remark 10 If the tensor is given by f l60|) . then the (initial) directors belong to 
the tangent plane at any point of and d^ = n°. Indeed, let P = Q^U be the polar 



decomposition of P. Using the matrices of components in the {cj ® ej} tensor basis, 
we write this relation as P = Q^U^ , and from (140 p we derive consecutively 
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where u^^p andvP^p are some given real functions of {xi,X2)- In view of fllip . it follows 

^ g°e3 = nO, (61) 

from which we can see that the third column of the matrix is equal to nP . On the 
other hand, by the definition (^2, the initial rotation field is given by = d^^Si 
and the matrix can be written in column form as 

4) . (62) 

/ 3x3 

If we compare ( 1611) and ( |62l) we find that d'^ = . Thus, we have = 'nP and {<i°, dg} 
is an orthonormal basis in the tangent plane, at any point of S^. 

If we choose the tensor as in fl60|) . then in order to satisfy ([33]) 2 we need to 
consider an additional regularity assumption on the initial configuration, namely 

polar(P) = polar (V0(xi,X2,O)) G H\uj,S0{3)), 

which is equivalent to Curl [polar(V0(xi, 0:2, 0))] G L^{u, S0{3)) , cf JJ^ . A stronger 
sufficient condition is & e W^'°°{uj,R^) n H'^{uj,R^). 

It is possible to simplify the form of the equations in the case of an orthogonal 
parametrization of the initial surface 5**^ . If we assume that the curvilinear coordinates 
{xi, X2) are such that the basis {ai, 02, n°} is orthonormal, then the initial surface 
is formally parametrized by orthogonal arc-length coordinates [13] and we have 

tta = a", aa(3 = a°'^ = ^ap ■ (63) 



Q° = PU' 



0,-1 



ai 



a2 



n 



u\i 



U\2 



u 



3x3 



12 



22 












Remark 11 We mention that such orthogonal parametrization satisfying the relations 
fl63|) exists only for developable surfaces 5° , since the Gaussian curvature is necessarily 
zero in this case. For example, let be a cylindrical surface with generators parallel 
to 63 . The position vector is given by 

9 9 

= y^(9, z) = Tq cos — ei + ro sin — e2 + z (rg > constant). 



Choosing the curvilinear coordinates Xi = 6,X2 = z, we have 

r^n.^ & nn n ^ ■ ^ 

o-i = 0*1?/ = — sm — ei + cos — 62, a2 = = £3, n = cos — Ci + sm — 62 , 

ro ro ro ro 

so i/iai {01,02,71"} is orthonormal. 

In view of (I60i) -( l63l) . we obtain in this case that = P (since C/° = II3 and 
polar(P) = P e S0{3)) and the directors {d^} in the reference configuration coincide 
with {ai, a2, n°} in each point of , i.e. 

dl = a^ = d^y\ d0 = a3 = n°, = Q V , d, = RQ^ ^ai. (64) 

The expressions of the elastic strain measures E'^ and may be simphfied in this 
situation. By virtue of (1^ and we get 

Q° - Q^'^d^y') = Q'^R^d^y - a„ ^^^^ 

= [ai ■ Q^R^d^y - 6ai)a, = [di ■ d^^y - 6^i)a, . 

We can write R^daR = {di ® ei)'^{dadj (g) ej) = {di ■ dadj)ei ® ej , so that we find 
axl{R^daR) = ^ eijk{dk ■ dadj)ei , axl(Q°'^(9aQ°) = ^ eijA:(afc ■ daaj)ei , 

where Cj^fc is the permutation symbol. The last relations and Q^''^ai = yield 

Q° raxl(il^9„il) - axl(Q°'^9^Q°)l = r(axl(i2^9«i?) - axl(Q°'^9„Q°)) -6^1 

1 r 1 (66) 

= 2 [{^k ■ dadj) - {ak ■ <9a%)J ai . 

Using the relations (IT7|) . ( 12^ . (1^ and ( 1^ we decompose the strain tensor and 
the curvature tensor in the basis {ai® a^} as follows 

= Ef^ai ® a^ , Ef^ = ■ (9„?/ - , = i^^a^ ® , 

kf^ = d3-dad2-a3-daa2, kl^ = di-dad-i-ai-daas, K^^ = d2-dadi-a2-daai. 

(67) 

For later reference, we introduce the notations 

El = E^ - {n^ ®n^)E^, K\ = - {n^^ ® 'nP)K\ (68) 
Then, from ( 167|) we get 

= El^a^ ® , no^;^ = E^^^n^ = E|„a„ , 



If we denote the matrices by = (EfA„ „ , K'^ = (KfA^ , , and also 
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(71) 



These expressions are completely similar to the strain measures for planar-shells in- 
troduced in pg|li^. 

Let us discuss next some important classes of elastic shells. 



5.1 Isotropic shells 

In the resultant 6-parameter theory of shells, the strain energy density for isotropic 
shells has been presented in various forms. The simplest expression of W{E'^, K^) has 
been proposed in the papers p^ [T5] in the form 

2W{E', K')= C[u (tr£;p2 + {i - y) ti{E''^^ E\)\ + a, C(l - u) n'^E^E^'^n^ 

+ D [u (trKp2 + (1 - z/) tr(K]j'^K|j)] +atD{l-iy) n'^K^K^^'^n^ 

(72) 

where C = is the stretching (in-plane) stiffness of the shell, D = i2{i^iy2) 
bending stiffness, h is the thickness of the shell, and , at are two shear correction 
factors. Also, E and u denote the Young modulus and Poisson ratio of the isotropic 
and homogeneous material. By the numerical treatment of non-linear shell problems, 
the values of the shear correction factors have been set to = 5/6, at = 7/10 in 



The value = 5/6 is a classical suggestion, which has been previously deduced 
analytically by Reissner in the case of plates [2S1IS1]- Also, the value at = 7/10 was 
proposed earlier in [501 see p. 78] and has been suggested in the work [19]. However, the 
discussion concerning the possible values of shear correction factors for shells is long 
and controversial in the literature [55115^] . 

With the help of the matrices (1701) . we can express the strain energy density (I72p 
in the alternative form 



e,T^0||2 



2W{E^,K'') = 

= C(l-z/)(||dev2symE||f + ||skew f ) + C^(tr \\E 

+D(l-z/)(||dev2symir,| f + ||skewir|| f ) + D^-^ [ti + at D{l-u) 

^ (73) 
where sym X = | (X + X"^) is the symmetric part, skew X = ^ (X — X"^) is the skew- 
symmetric part, and dev2X = X — i(trX)lL2 is the deviatoric part of any 2x2 
matrix X. The coefficients in (173!) are expressed in terms of the Lame constants of the 
material A and fi by the relations 

2 2/i + A ' ^ ' ^ ' 2 12 2/i + A ' ^ ' 6 

Then, we obtain that the given quadratic form (173|l is positive definite if and only if 
the coefficients E and u satisfy the inequalities 

E>0, -1 < z/ < ^ . (74) 

In terms of the Lame moduli of the material, the inequalities ([71]) are equivalent to 

/x>0, 2/i + 3A>0. 

These conditions are guaranteed by the positive definiteness of the 3D quadratic elastic 
strain energy for isotropic materials. Thus, we find that the strain energy W is convex 
and satisfies the coercivity condition f l35|) . so that the hypotheses of Theorem |6] are 
fulfilled. Applying Theorem [HI we obtain (under suitable assumptions on the given 
load and boundary data, and the reference configuration {y'^,Q^)) the existence of 
minimizers for isotropic shells with strain energy density in the form (!72|) . 

In [21], Eremeyev and Pietraszkiewicz have proposed a more general form of the 
strain energy density, namely 

2W{E',K') =ai{tTElY + a2tT{ElY + a3tr{El'^El) +ain'^E'E''^n'^ 

+/3i(trK|jf + f32tT{KlY + /33tr(Kj'^Kf) + fS^n^K^K^'^n^ . ^ ' 

The eight coefficients , /3k {k = 1, 2, 3, 4) can depend in general on the structure cur- 
vature tensor = ax\(^daQ^Q^'^) of the reference configuration. For the sake of 



simplicity, we assume in our discussion that the coefficients and /3k are constant. We 
can decompose the strain energy density (175]) in the in-plane part Wpig^ne{E^) and the 
curvature part Wcnvv{K'^) and write their expressions using the matrices of components 
rOl) in form 

WiE', K') = Wp,,UE') + W,^,,{K') , (76) 



2Wpiane(£^') = (fts + c^s) || sym Ef, f + (ftg - i| skcw K f + «i (trE,^) ' + «4 1| ^ 



2W,^r,iK') = (/32 + /33)||symirff +(/33-/32)||skewK|f +/3i(trir|)' + /34||K 
The in-plane part of the energy density (1761) can be written equivalently as 

2H^piane(-E^) = (^2 + c^s) || devs sym E| \\^ + {a^ - ^2) || skewEf f 

V ' V ' 

in-plane shear-stretch energy in-plane drill rotation energy 

+ (a, + ^^){tr^Y + a,\\E-^n^r . ^^^^ 

^ transverse shear energy 

in-plane elongational stretch energy 

The above forms of the strain energy W are expressed in terms of the components of 
the tensors and K'^ in the basis {a^ 0.0} , i.e. in terms of the elements of the 
matrices ( ITOl) . Denoting with /i^"" ^^e coefficient (0:3 — 02) in ( ITTl) - we remark that 
the term 

/if" II skew ^1 f , with /if" := as - , (78) 

describes the quadratic in-plane drill rotation energy of the shell. We call the coefficient 
^driu linear in-plane rotational couple modulus, in analogy to the Cosserat couple 
modulus in the three-dimensional Cosserat theory [39]. 

Remark 12 The planar isotropic Cosserat shells have been investigated also in 138^4'^ ! 
using a model derived directly from the 3D equations of Cosserat elasticity. We mention 
that the expressions ( 1761) . ( 1771) of the strain energy density are essentially the same as 
the strain energy of the Cosserat model for planar shells 138]. By comparing these two 
approaches (6-parameter resultant shells and Cosserat model) we deduce the following 
identification of the constitutive coefficients ai , a4 

ai = h "^^^ , a2 = h{fi- He), cts = /i(/i + /i^), ^4 = k /i(/i + /ij, (79) 
z/i + A 

where fic is the Cosserat couple modulus of the 3D continuum, and k is a formal shear 
correction factor. From ( 175]) . (17^ we observe that 

/if" = a3-a2 = 2/i/i,, (80) 



^n0||2. 



which means that the in-plane rotational couple modulus /if ^/ Cosserat shell 
model is determined by the Cosserat couple modulus fic of the 3D Cosserat material. 



The relations fl79p are similar to the corresponding relations in the linear theory 
of micropolar plates, see 13, Eqs.(45)]. From a mathematical viewpoint, the difference 
between the two sets of relations consists in the notations used and the value of the 
shear correction factor. 

Looking at f lTB]) and f l77|) we observe that the quadratic form W{E^, K^) is positive 
definite if and only if the coefficients verify the conditions 

2ai + ^2 + as > 0, 0:2 + as > 0, as - a2 > 0, a4 > 0, 
2/3i + /32 + /3s > 0, /32 + /3s>0, /3s - /32 > 0, /34 > 0, 

Provided that the conditions (1811) are satisfied, the strain energy function W{E^, K^) 
is convex and coercive in the sense of fl35|) . By virtue of Theorem [6l in this case the 
minimization problem associated to the deformation of isotropic elastic shells admits 
at least one solution. 

Remark 13 The same conditions (IHTj) have been imposed in ^22] to establish existence 
results in the linearized theory of micropolar (6-parameter) shells. 

Remark 14 The case /i^"^^ = (i.e., — a2 = is not uniformly positive definite. 
However, with a slight change of the resultant shell model, one can prove the existence 
of minimizers using similar methods as in j^Bj. A linearization of such a model leads 
exactly to the Reissner kinematics with 5 degrees of freedom fj^ , where the in-plane 
drill rotation is absent. The physical meaning of the in-plane rotational stiffness /x^"" = 
as — a2 in the resultant shell model is not entirely clear to us. 

Since only two independent rotations are required to orient a unit director field, a 
distinctive feature of classical plate and shell theories is a rotation field defined in terms 
of only two independent degrees of freedom. Rotations about the director itself - the 
so-called drill rotation - are irrelevant and for that matter undefined in classical shell 
theory. 



5.2 Orthotropic shells 

The constitutive equations for orthotropic shells have been presented in [21] within 
the 6-parameter resultant shell theory. The expression of the strain energy density in 
terms of the tensor components defined in fl67I) is given by 

2W{E\ K^) = EI^E;, + Df^ El^% + C^p^, K^pK^s + D% KI^K^^ (82) 

where C^^^^ , C^^^ , D^^ and are material constants which satisfy the following 
symmetry relations 



We observe that the quadratic function (1821) is coercive if and only if the following 
symmetric matrices are positive definite 



111 

'-^1122 
'-^1112 



1121 



riE 
'-^1122 

/-iE 
'-"2222 

r^E 

nE 

0< 



2221 



riE 
'-"1112 

r-iE 
'-'2212 

r^E 
'-"1212 
r'E 



1221 



/^E ■ 

'-"1121 

^E 
'-"2221 

r^E 
'-"1221 
r'E 



2121J 



•-"1111 


'-"1122 


'-"1112 


c 


'-"1122 


r<K 

'-"2222 


rK 

'-"2212 


Q 


'-"1112 


'-"2212 


rK 

'-'1212 




r<K 

.'-"1121 


r'K 

'-"2221 


rK 

'-'1221 





K ■ 
1121 

K 

2221 
1221 
2121. 



^12 



12 

r)E 

^22 



jjK jjK 

-^11 -^^12 

r)K jjK 

'^Vl ^22 



(83) 

In the situation when the matrices (1831) are positive definite, then the strain energy W 
given by f l82|) satisfies the hypotheses of Theorem [61 Then, we can use our theoretical 
results to derive the existence of minimizers for orthotropic shells. 



5.3 Composite layered shells 



Let us analyze the case of composite shells made of a finite number of individually 
homogeneous layers. According to [13j, the strain energy density of such type of shells 
can be written by means of the tensor components fl671) in the form 



4) 



where A^p^s , Ba/s-ys , Dap-ys , Sajs and Gap are the constitutive coefficients of composite 
elastic shells, which have been determined in |13] in terms of the material/geometrical 
parameters of the layers. They satisfy the symmetry conditions 

^a(3^S = ^7<5a/3 ; D^fi^S = D^^^P i Sap = S^a , Ga/3 = G pa ■ 

In the constitutive relation (184]) one can observe a multiplicative coupling of the strain 
tensor E'^ with the curvature tensor for composite shells. Let us denote by A, D 
and B the 4x4 matrices of material constants 



A 



Ann 


An22 


An\2 


^1121 






I^iiii 


-D1122 


Dni2 


-D1121 


An22 


A2222 


^2212 


^2221 




D = 


-D1122 


-D2222 


-D2212 


-D2221 




















^1112 


^2212 


^1212 


^1221 






-D1112 


-D2212 


-D1212 


-D1221 


^1121 


^2221 


^1221 


^2121. 






.-D1121 


-D2221 


-D1221 


-D2121 








-Bun 


-B1122 


Bni2 


-B1121 












B = 


-B2211 


-B2222 


-B2212 


-B2221 














-B1211 


-B1222 


-S1212 


-B1221 














.-B2111 


-B2122 


-S2112 


-B2121. 









One can show that the necessary and sufficient condition for the coercivity of the strain 
energy function (1841) is that the following matrices are positive definite 



C 



v44x4 -84x4 
-84x4 -^4x4 



8x8 



'S'll >S'i2 
'S'12 >S'22 



2x2 



Gn G12 
G12 G22 



2x2 



With these notations, one may write the strain energy density f l84p in the matrix form 

2W{E^, K^) = VCV^+{ El, , El^) § ( El, , El,f + ( Kl, , Kl,) G ( Kl, , Kl,)^ , 
with V = ( Ef, , E^,, , E^,, , E^,, , ir.^i , ^^2^2 , Kf, , • 

In conclusion, if the matrices C, § and G are positive definite, then we can apply 
Theorem [6] for the strain energy density given by fl8^ and prove the existence of 
minimizers for composite layered shells. 

Remark 15 The results and conclusions presented above are obviously valid also in 
the case of plates, i.e. when the reference base surface is planar. However, many of 
the formulas for general shells can be significantly simplified in the case of plates, since 
the 3 orthonormal bases {ai,a2,n°}, {d^,d2,d^} and {61,62,63} can be considered 
identical. 

The corresponding existence results for 6-parameter geometrically non-linear plates 
(planar shells) has been presented in 111] for isotropic or anisotropic materials, and 
in fW^ for composite planar-shells. We mention that, in the case of isotropic plates, 
the existence theorem can be obtained from the more general results concerning Cosserat 
planar-shells presented by the second author in f3^\4^. 

In a forthcoming contribution we will extend our existence results to the 6-parameter 
resultant shell model with physically non-linear behavior and show the invertibility of 
the reconstructed deformation gradient F. 
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